There are a number of reasons for organizations and governments to promote the reduction of meat consumption, including the goals of reduced animal suffering, reduced rate of climate change and fresh water contamination, food sustainability, and the potential improvements in population health. As a result, different campaigns aimed at this goal have been employed by various environmental, animal rights, and health focused groups. One such campaign, sometimes called “Meatless Monday”“, focuses on participants eliminating eating animals for one day per week. Some of these campaigns have made claims of how much water, greenhouse gas emissions or number of animals that this participation is expected to save; however, these estimates assume both that a participant’s meat consumption is consistent with average meat consumption and that a participant’s meat consumption on days other than Monday (or whatever day chosen) does not change due to their participation in the campaign.
It is not currently clear how reasonable these two assumptions are. It may be that the meat consumption of those open to participating in these programs is not consistent with the average level. These people are likely to have prior knowledge of some of the adverse effects of eating meat advertised by the campaign, so a person with an interest in participating in one of these programs may already be eating less meat than average. It also may be reasonably argued that actively choosing to reduce meat consumption on one day of the week could result in the participant modifying their behavior on the other days of the week. The participant could consciously or unconsciously compensate for their meatless day by eating more meat on other days of the week. Alternatively, participation on a weekly basis may keep the goal of reduced meat consumption more present in the participant’s mind and help them to further reduce meat consumption on other days of the week. This experiment aims to quantify the true effect on overall meat consumption, as this would be useful to organizations interested in promoting reduction in order to compare the value of such a campaign against other areas of investment.
The study follows a two-group, within-subjects design as depicted in Figure 1. The primary data collected was the weight of each participant’s meat intake over the course of the study. All subjects were observed for one week to gather daily measurements of meat consumption through self-reported surveys. In order to help ensure the measured effect is an unbiased estimate of the effect of treatment, subjects were then randomly divided into treatment and control groups. The night before the beginning of the second week of observation, the participants randomized into the treatment group received an additional email asking them refrain from eating meat one day in the upcoming week. They were asked to identify the day of their choice using an additional single question survey. No additional contact was made with participants in the control group.
Figure 1: Experimental Design
Participants were asked to weigh meat or examine meat packaged for the given weight, if possible, and failing this provide a description along with a weight or volumetric estimate. Accepting estimates in addition to measured data was expected to improve compliance relative to a more strict procedure. In order to aid participants and promote consistency among estimates, the same guidelines on judging meat weight were provided at the beginning of each daily survey. These guidelines are provided in Figure 2.
Figure 2: Example of Survey Guidelines
Potential study participants were obtained through an online Qualtrics survey. Links were sent out identifying that a study was recruiting participants interested in reducing the amount of meat they were currently eating. Only respondants that selected that they had an interest in this were accepted in the study. The effect of treatment is only desired for the subset of the population interested in reducing their meat consumption, since Meatless Monday campaigns are voluntary. The effect among the general population, if it could be obtained, would be of less value. This survey was also used to obtain demographic information on each participant that was thought could be potentially useful in identifying heterogeneous treatment effects.
As shown in Figure 3 outlining the experimental procedure, 54 responses were obtained from the initial survey. Fifty of these people identified that they were interested in reducing the amount of meat they were eating and signed up for the experiment. Of these, 12 did not respond to any further follow-up data collection emails.
Figure 3: Experiment Progression
Figure 4 shows the attrition and non-compliance. There are 29.5% attrition before the first week, and 18.4% attrition after the first week. Most of the attrition occured before the subjects were randomized into treatment and control groups, and can be definitively said is not due group assignment. Therefore attrition is not expected to introduce a systemic bias. In order to avoid the effect of attrition on randomization, retain subjects and attrite subjects are randomly assigned treatment or control separately. After the second week, 6.9% attrition happened. One subject missed reporting one day in the first week and one day in the second week.
There are 12.5% of those in treatment group did eat meat on the day they selected to abstain. This is likely to occur during in meatless campaigns outside experimental conditions as well, and does not present a problem. The results from those in a meatless campaign who choose to eat meat are of interest as well. The primary measure of interest is the “Intent to Treat” of a meatless campaign among those interested in order to examine the true effectiveness of these programs.
Figure 4: Bar chart and table hierarchy chart shows the aggerated attrition and non-compliance number. The table hierarchy also noted the experiment design and randomization process. Analysis all based on the orange data set. .
The data provided in Figure 5 shows that subjects in control have more variability in daily meat consumption than those in the treatment group from a daily meat consumption perspective. Across both the weeks, the average daily meat consumption for the treatment group is less than that for the control group by about 2.4 ounces.
Figure 5: Distribution of daily meat consumption by treatment and control groups
Figure 6 provides the distribution of meat consumption for 2 weeks for both the treatment and the control group. Both the treatment and control group are in “control” for week 0 and the treatment group goes in to treatment in week 1. For both weeks, the average daily meat consumption per person for the control group is greater than that of treatment group. Going in to week1, the mean for the control group has increased by about 0.4 ounce and the mean for treatment group has decreased by about 0.4 ounce. This suggest that from a baseline meat consumption perspective, people in the treatment group already consume less meat than those in the control group.
Figure 6: Distribution of daily meat consumption by treatment and control groups
Figure 7 provides meat consumption variation over different stages of treatment. Average meat consumption per day for the treatment group under control conditions during the first week is 5.56 ounces. Once the treatment group is asked to refrain from eating meat on a day of their choice, there may be some anticipatory effect. During the time the treatment subjects are waiting in anticipation of the treatment day, their average daily meat consumption goes up by near 1 ounce. On the day of the treatment, a reduction in meat consumption is expected due to the nature of treatment. The average daily meat consumption reduces to 2.25 ounces for those in treatment on the day. After the treatment day, their average daily meat consumption increases to 5 but does remain lower than the baseline average. However, it should be noted that some subjects elected the first day of the week and some chose the last. Given the short window of treatment, pre-treatment and post-treatment data is not available for all those in treatment. Additionally, the number of days in the pre and post treatment day periods vary among subjects. Therefore, care should be taken to avoid drawing conclusions from these comparisons.
Figure 7: Distribution of daily meat consumption for the treatment group by treatment stages
Figure 8 provides the daily responses by subject. There are three subjects that are driving up the average daily meat consumption in the control group. The subjects have high variability as far as their daily meat consumption goes. For the treatment group, there is one subject driving up the average daily meat consumption.
Figure 8: Individual Distribution of daily meat consumption
The initial survey gathered information in eight different categories. The ‘Interest_Level’ covariate is to check the subject’s interest level of reducing meat consumption. Study requirements removed the 4 subjects not interested in reducing meat consumption, leaving subjects who selected either moderately interested or very interested. A covariate labeled ‘Reasons’ asked people their primary reason they were interested in to reducing the amount of meat they eat. It is possible that people different motivations cause people to react differently to treatment. ‘Household_Vegetarian’ variable is to check if there is any vegetarian in household. People who have a vegetarian family member may have different reaction to the treatment. ‘People_in_Household’ covariate is to present the number of people in household. Since people in a single household tend to eat meals togther, those with larger families may have less influence on daily meals due to other family members cooking or shopping. ‘Sex’, ‘Age’, ‘Income’, and ‘Education’ are the demographic covariates collected.
Figure 9 shows the distribution of the covariates in treatment and control group. Differences between treatment and control are likely the result of the low number of subjects in some of the groups rather than an issue with the randomization process.
Figure 9: Summary of proportion and bar chart of covariates
Figure 10 provides the results of the t-test and chi-squared test examining covariate balance between the treatment and control group. The null hypothesis is that there is no difference in proportions of the covariates from treatment and control group. The p-value of all the covariates are above 0.05, indicating that we cannot reject the null hypothesis. Notice that the covariate ‘Reasons’ has the least p-value 0.1332, which is mainly due to the unbalanced ‘Reduce Animal Suffering’ mentioned previously.
Figure 10: Covariates Balance Check by t-test and chi-squared test. P-values > 0.05 means that no information indicating failed randomization
A simple regression was performed on the data from week one only during which all subjects were in control to examine whether a baseline difference in meat consumption was present between those assigned to treatment and those assigned to control in the second week. The dependent variable here is the daily meat weight measurement and the independent variable is an indicator variable for whether the observation comes from a subject assigned to treatment in the second week of the experiment. Although the difference is not statistically significant, it shows our best estimate is that on average those assigned to the treatment group eat approximately 2 ounces less meat per day than the subjects that were assigned to the control group for the duration of the experiment. It should be noted that all estimated standard errors and measures of uncertainty reported within this section are cluster-robust to account for the fact that each of the daily observations from a subject are related to each other.
##
## Table 1: Simple Regression for Baseline Differences
## ======================================================
## Dependent variable:
## ----------------------------------
## Daily Weight of Meat Consumed [oz]
## ------------------------------------------------------
## Treatment_Group -2.082
## Constant 7.644*** (5.012, 10.276)
## ------------------------------------------------------
## Observations 202
## R2 0.021
## Adjusted R2 0.016
## Residual Std. Error 7.049 (df = 200)
## F Statistic 4.353** (df = 1; 200)
## ======================================================
## Note: *p<0.1; **p<0.05; ***p<0.01
Although not a statistically significant difference, this indicates it will be important to control for these baseline differences in tendency when estimating the causal effect of treatment overall. The next regression shows that we do see a large, statistically significant effect on the amount of meat eaten on the treatment day due to treatment assignment. This is consistent with expectation since most of those assigned to treatment did follow through and report zero meat consumption for that day.
##
## Table 2: Regression of Meat Consumption on Treatment Day Indicator
## ======================================================
## Dependent variable:
## ----------------------------------
## Daily Weight of Meat Consumed [oz]
## ------------------------------------------------------
## TreatmentDay -3.356
## factor(ID)2 -12.714
## factor(ID)3 -10.500
## factor(ID)4 -0.429
## factor(ID)5 -5.429
## factor(ID)6 -12.357
## factor(ID)7 -6.286
## factor(ID)8 -8.929
## factor(ID)9 -8.857
## factor(ID)10 -6.143
## factor(ID)11 -6.429
## factor(ID)12 -7.571
## factor(ID)13 -13.679
## factor(ID)14 -14.689
## factor(ID)15 -6.117
## factor(ID)16 -8.546
## factor(ID)17 -11.403
## factor(ID)18 -11.546
## factor(ID)19 -9.332
## factor(ID)20 -14.903
## factor(ID)21 -12.832
## factor(ID)22 -13.546
## factor(ID)23 -7.617
## factor(ID)24 -13.975
## factor(ID)25 0.668
## factor(ID)26 -4.403
## factor(ID)27 -7.689
## factor(ID)28 -11.332
## factor(ID)29 -9.903
## Constant 15.429*** (15.429, 15.429)
## ------------------------------------------------------
## Observations 404
## R2 0.430
## Adjusted R2 0.386
## Residual Std. Error 5.124 (df = 374)
## F Statistic 9.731*** (df = 29; 374)
## ======================================================
## Note: *p<0.1; **p<0.05; ***p<0.01
However, this does not answer the primary question of whether a difference is maintained over the course of the week. The next regression examines the treatment effect where a treatment observation is defined as one that occurs within the second week from a subject that was asked to refrain from eating meat one day that week. This regression indicates that the best estimate when controlling for individual behavior is that those subjects reduce overall meat consumption by an average of 0.39 ounces each day of the week they are in treatment. However, the 95% confidence interval on this effect is quite large, indicating that the experiment does not have the statistical power to determine whether an effect of this size is simply due to random variation. This indicates that although subjects in treatment tended to comply with not eating meat on the day of their choice, they balanced this with eating more meat than they normally would on other days in that week. Although the best estimate of the difference is that those in treatment maintained a small reduction, more subjects would be needed to reduce the uncertainty and confirm this.
model_treatment <- lm(Weight~InTreatment+factor(ID),data=all)
stargazer(model_treatment, ci=TRUE, se=sqrt(diag(cluster.vcov(model_treatment,(factor(all$ID))))),header = FALSE, no.space = TRUE,single.row = TRUE,type='text',title = "Table 3: Regression of Meat Consumption on Treatment",
dep.var.labels=c("Daily Weight of Meat Consumed [oz]"))
##
## Table 3: Regression of Meat Consumption on Treatment
## ======================================================
## Dependent variable:
## ----------------------------------
## Daily Weight of Meat Consumed [oz]
## ------------------------------------------------------
## InTreatment -0.393
## factor(ID)2 -12.714
## factor(ID)3 -10.500
## factor(ID)4 -0.429
## factor(ID)5 -5.429
## factor(ID)6 -12.357
## factor(ID)7 -6.286
## factor(ID)8 -8.929
## factor(ID)9 -8.857
## factor(ID)10 -6.143
## factor(ID)11 -6.429
## factor(ID)12 -7.571
## factor(ID)13 -13.679
## factor(ID)14 -14.732
## factor(ID)15 -6.161
## factor(ID)16 -8.589
## factor(ID)17 -11.446
## factor(ID)18 -11.589
## factor(ID)19 -9.375
## factor(ID)20 -14.946
## factor(ID)21 -12.875
## factor(ID)22 -13.589
## factor(ID)23 -7.661
## factor(ID)24 -14.018
## factor(ID)25 0.625
## factor(ID)26 -4.446
## factor(ID)27 -7.732
## factor(ID)28 -11.375
## factor(ID)29 -9.946
## Constant 15.429*** (15.429, 15.429)
## ------------------------------------------------------
## Observations 404
## R2 0.421
## Adjusted R2 0.376
## Residual Std. Error 5.165 (df = 374)
## F Statistic 9.372*** (df = 29; 374)
## ======================================================
## Note: *p<0.1; **p<0.05; ***p<0.01
#coeftest(model_treatmentday,vcov. = cluster.vcov(model_treatmentday,all$ID))
As noted previously, the behavior while under treatment may vary on when in the week the subject is relative to treatment, but due to the limited duration in treatment, data for a consistent comparison across the treatment group is not available. This kind of data would be valuable in future studies.